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I. INTRODUCTION 



Perturbations of black holes are dominated, over intermediate timescales, by characteristic modes known as quasi- 
normal modes. They are similar to normal modes of a closed system, but since the field can fall into the black hole 
or radiate to infinity, the modes decay and the corresponding frequencies are complex 0. They have been exten- 
sively studied in asymptotically flat spacetimes and good reviews on the topic can be found in 2], 3]. The study of 
quasinormal modes in Anti-de Sitter space was first done in 0. 

More recently quasinormal frequencies have been investigated in the context of string theory in Anti-de Sitter space 
[l|i 0- There is a suggestion, known as AdS/CFT correspondence, that string theory in Anti-de Sitter space is 
equivalent to a conformal field theory in one dimension fewer 0. In the framework of this conjecture the study of 
AdS black holes has a direct interpretation in terms of the dual conformal field theory on its boundary. The duality 
predicts that the retarded CFT correlation functions are in one to one correspondence with Green's functions on 
Anti-de Sitter space with appropriate boundary conditions 0, 01 > E3> E2 Furthermore, as mentioned in 0, it 
is assumed that a large static black hole in AdS space corresponds to a thermal state in the CFT on the boundary. 
Perturbing the black hole is equivalent to perturbing this thermal state. The perturbed system is expected, at late 
times, to approach equilibrium exponentially with a characteristic time-scale. This time-scale is inversely proportional 
to the imaginary part of the poles of the correlators of the perturbation operator. It seems that these relaxation time- 
scales are quite complicated to calculate in the CFT, therefore their computation is conveniently replaced by the 
evaluation of the quasinormal frequencies in the AdS bulk space. In and 0, the authors went beyond the scalar 
perturbation treated in and they considered electromagnetic and gravitational perturbations. In one of the most 
recent works on the subject, Berti and Kokkotas 13] confirm and extend previous results on scalar, electromagnetic 
and gravitational perturbations of static AdS black holes and analyse, for the first time, Reissner-Nordstrom Anti-de 
Sitter black holes and calculate their quasinormal frequencies. 

The interaction of a Dirac field with a black hole has been studied by Finster and his collaborators in a series of 
papers. While they have found stable particle- like solutions in the Einstein-Dirac- Maxwell system [Til Il5|] . they have 
also proved the non-existence of time-periodic solutions in various black hole space-times |l6l Il7|| . This means that 
Dirac particles, including electrons and neutrinos, cannot remain on a periodic orbit around a black hole. In [Tsj]. they 
also showed that there are no spherical symmetric black hole solutions in the Einstein-Dirac-Maxwell system other 
than the Reissner-Nordstrom one. This suggests that if a cloud of Dirac particles undergoes gravitational collapse, 
the fermionic particles either vanish inside the event horizon of a black hole or escape to infinity. 
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The quasinormal frequencies related to the evolution of a massless Dirac field in a Schwarzschild black hole space- 
time, were studied in [lij. The quasinormal modes of the Reissner- Nordstrom de Sitter black hole for Dirac fields 
were studied by using the Poshl- Teller potential approximation in Ref. j2lj. In this paper we are interested in the 
analysis of the modes of vibration of a massless Dirac field in a Schwarzschild Anti-de Sitter background space and 
it is for this reason that we compute the related quasinormal mode frequencies. This calculation is relevant to the 
AdS/CFT correspondence between M-theory on AcLSa X S 7 and SU(N) super Yang-Mills theory on M3, and it also 
serves to fill in a gap in the literature, which, in AdS spaces, has only been concerned with particles of integral spin 
0,1,2. 



II. DIRAC EQUATION IN STATIC ADS SPACE-TIME 



The metric for a Schwarzschild Anti-de Sitter black hole can be written as 



ds 2 = -fdt 2 + jdr 2 + r 2 {d0 2 + sin 2 edcp 2 ), (2.1) 

with 

/ = l-^-^, (2.2) 
r 3 

where the parameters M, and A represent the black hole mass, and the negative cosmological constant, respectively. 
The Dirac equation in a general background space-time can be written, according to |2l| , as 

[7 a e a na M + r A1 )]* = o. (2.3) 

Here, 7° are the Dirac matrices, 
while a % are the Pauli matrices, 



1 f0 1 \ 2 / -i\ 3 / 1 
° = ' 1 ' a = [i ' ° = -1 



(2.5) 



The four- vectors e a M represent the inverse of the tetrad e a M defined by the metric as, 

9tiv = ?7abe a M e b l/ , (2.6) 
with r/at) = diag(— 1, 1, 1, 1) being the Minkowski metric. are the spin connection coefficients, which are given by 

r M = i[7 a ,7 6 ]e "e 6 „ il ,. (2.7) 

Here, = d^e^v — I 1 " ei, a is the covariant derivative of and r" is the Christoffel symbol. 
We take the tetrad to be 

e£ = diag(V7, -^,r,r sin <?). (2.8) 
The spin connection can then be expressed as 

To = j f loll, 

Ti = 0, 

r 2 = 2^/7172, 

T 3 = i(sin 6»V/7i 73 + cos ^7273). (2.9) 
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The Dirac equations l|2.3|) become 



70S* r— f d 1 1 df\ T 72 ,9 1 „ T 73 9* „ 
V / ^ \or r 4j dr J r 06 2 r sin a<£ 



(2.10) 



If we re-scale * as 

* = /-3$, (2.11) 
Eq. (|2.1UI) assumes a simpler form in the new unknown <I>, which can be written as 

-^f + ^(! + ;)* + t<! + 5 -•>• + ^ - <"*> 

We introduce a well known coordinate change from the radial variable r to the tortoise coordinate r* given by 



We will use an ansatz for the Dirac spinor 



4>f m (0,V>) 



*= I P </>^T\ ' l e > ( 2 - 14 ) 



with spinor angular harmonics 



— ' ,„ I , {for 3 = 1+ -), (2.15) 



2j"+2 J i 

Note that K m±1/ ' 2 (#, (p) represent ordinary spherical harmonics. Since 



j+i-m ym-1/2 \ 

2j+2 ' , , (/or 7 =i--). (2.16) 



-i(& + |cat 

sin 6 c?ip 




(2.17) 



ea uat ions 12 . 1 21 can be written in the simplified matrix form 



o-„WF^ a ^V^f* )(£]=„. ("si 



■ ;v G v"*:VG ± y i TVJ i, -^A G 

These equations can be thought of as equations referring to the radial functions (F + , G + ) or (F~ , G~) depending on 
the original choice made in the ansatz of eauation l2.14l On the other hand it can be proved that the radial differential 
equations related to the two different sets of functions are exactly the same. In other words, we have a unique equation 
for i r± and another equation for G^. For this reason we will avoid from now on specifying if we are referring to the 
spin up or spin down components when talking about and and we will simply call these functions F and G. 
Consequently the two decoupled equations can be expressed in the form 

^. + ( ( J a -V 1 )F = 0, (2.19) 

+ (uj 2 - V 2 )G = 0, (2.20) 
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with 



(*— i-i. >— i). 

In a similar way to the integer spin case in [Tsj . the two potentials V\ and V2, are super-symmetric partners derived 
from the same super-potential. In the asymptotically flat case, it has been shown that potentials related in this way 
possess the same spectra of quasinormal mode frequencies 22] . In the AdS case, this is still true provided that the 
boundary conditions are mapped between the two cases in a consistent way. We shall concentrate just on Eq. 12.191 
with potential V\ in evaluating the quasinormal mode frequencies in the next section. 

III. DIRAC QUASINORMAL FREQUENCIES 

In this section we evaluate the quasinormal frequencies for the massless Dirac field using the Frobenius series 
solution method. We start defining A as 

A = r 2 - 2Mr + a 2 r i , (3.1) 

where a 2 = —A/3. This implies that the function / in 12. 191 can be re-written as / = A/r 2 . Next, rescale F as 

F = e- iujr *u, (3.2) 

which allows us to scale out the behaviour of the function at the black hole event horizon. We obtain an equation in 
u given by 

„ 9 d 2 u ( df 1 du , 



dr 2 \ dr J dr 

Let us now set 

r = -. (3.4) 

x 

We obtain the equation in the variable x 

/V^ + fx 2 (2xf + x 2 ^- + 2iJ) ^-Vu = 0. (3.5) 
dx \ dx J dx 

Define p as 

p = x 2 - 2Mx 3 + a 2 x i . (3.6) 
Note that p\ — p{x\) = A{x\)x\ = 0. Rescale / as 

f = Ax 2 =px~ 2 (3.7) 

inserting this into 13.51 we obtain the equation 

2 d 2 u f dp n . \ du 



p d^ +p [Tx +2i " Tx- Vu = - (3 - 8) 



Rescale again the independent variable x, by introducing the variable z given by 

^^-)=z 2 , (3.9) 

Xl J 

we obtain a new expression for V in the variable z, 

V(z) = P \k\ 2 + P 1/2 \k\ Xl (1 - z 2 ) [3Afai (1 - z 2 ) - l] . (3.10) 
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Insert |3~§1 into f5~51 and rescale all the parameters: M, a 2 ,uj with respect to ii, which is the inverse of the radial 
coordinate r\ of the black hole event horizon, while r*i is the largest root of A. We obtain 

Pz ^ _ L _ g z + 4lt ^ 2 J ^ - 4z 3 {|fc| 2 + p- 1 / 2 [3M|fe| (1 - z 2 f - \k\ (1 - z 2 )] } u = 0, (3.11) 

where 

p = px^ 2 . (3.12) 

Expanding the coefficients of eauation l3.11l as Taylor series around z = 0, which is the same than performing a series 
expansion around the black hole event horizon, we obtain 

A(z)^-B(z)^-C(z)u = 0, (3.13) 
az^ az 



where 



with 



and 



A{z) = p 2 z 3 +p 3 z 5 + 2Mz 7 , (3.14) 
B(z) = {4iu-p 2 )z 2 -3p 3 z 4 -10 Adz 6 , (3.15) 
C(z) = 4\k\ 2 z 3 + 2\k\p 2 z 2 q- 1/2 +4\k\p 3 z 4 q- 1/2 + l2M\k\z 6 q~ 1/2 , (3.16) 



pa = (6M-2), (3.17) 
p 3 = (-6M + 1) (3.18) 

q=pz- 2 =p 2 +p 3 z 2 + 2Mz A . (3.19) 



The differential equation in it can be now solved by series by exploiting Frobenius method, which consists in looking 
for a solution of the form 



= z p J2 a nZ n - (3.20) 



u 

n=0 



The main difference between the spin ^ case and the integer spin case is the non-polynomial term q x l 2 . This function 
can be Taylor expanded as 



oo 

1/2 = (pa + P3Z 2 + 2Mz i )- 1 ' 2 = J2 ^sz s , (3.21) 



leaving aside for the moment the problem of deriving the coefficients of this last series. 

The two possible values for the index p at z = of our series solution 13.201 can be easily calculated and they are 
given by p\ = and p 2 — Aiuj/p 2 . We pick the first of the indices (the other solution has infinitely many oscillations 
close to the horizon) and therefore look for a series solution which will be simply given by 

oo 

u = a nz n (3-22) 

n=0 

The infinite nature of the expansion for q -1 / 2 and u, leads us to deal with a product of two series which can be 
expressed as a single series, 

(oo-. \ / oo \ oo / n \ 

E -1 A * ZS E = E E ~^n-s Z n (3.23) 

s=0 ' / \n=0 ) n=0 \s=0 ' / 
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Inserting equations 13 . 221 and 13 . 231 into eauation l3.13l we obtain a recurrence relation for a n+ i which has n terms, 
{(n + 1) [p 2 (n + 1) - 4iuj}} a n+1 + { P3 (n - 1) (n + 1) - 4|fc| 2 } a„_ x 
+ |2Af (n - 3) (n + 1)} a„_ 3 - 2\k\p 2 j A a„ + • • • + ^a | (3.24) 

- 4|fc|p 3 |A a„_ 2 + • • • + ^"^, «o| - 12M|fc| |A a„_4 + ■ ■ ■ + ^r^T a °} = 

This recurrence equation, due to its unbounded character caused by the series expansion for q^ 1 / 2 , is much harder 
to solve by numerical recursion, and in any case before we proceed we first need to know the expression for the 
coefficients A s of the Taylor expansion of q^ 1 ^ 2 . This problem is solved by remembering Faa di Bruno's formula for 
the analytical calculation of derivatives of any order 



d n h 



_dx n 

where h{x) = f[g(x)] is a composed function, 



= £/*E nI 11(0^1)^ (3 - 25) 



fk = —f(Vo) (3.26) 



dy 
and 

9i = -j-i9{no)- (3-27) 
The second sum inside l3~2*51 is done over partitions p(n, fc), which are defined as 

{n n ~\ 

(oi, . . . , On) : a, e No, £ a 4 = fc, £ ia* = n > (3.28) 
i=l t=i J 

where Ao is equal to the set of non-negative integers. An element (a±, . . . ,a n ) belonging to p(n,k) represents a 
partition of a set with n elements into a\ classes of cardinality 1,. . ., a n classes of cardinality n. The number of such 
partitions is represented by the Stirling numbers of the second kind, 

n , 

= E "'n^u-Twr-- (3-29) 
p(n,k) i=l 1 4 - ){ ■' 

When formula 13.251 is applied to the A s coefficients, we have 

A s = ± »(0)] £ Slfl—^^HO)}" , (3-30) 

m=l p(s,m) i=l V *' A ' y 

with / (m) [g(0)] given by 

"ff (-^f 11 ) ^(0)]- 1/2 - m - (3-31) 

On the other hand, evaluating q and its derivatives at z — 0, we obtain 

g(0) = 6M-2, (3.32) 

g"(0) = 2(1 -6M), (3.33) 

q( 4 )(0) = 48Af (3.34) 

while 

g'(0) = q"'(0) = 0, (3.35) 
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beside all the derivatives of order higher than four. The A s coefficients can then be written in the form 

A ° = £ /(m) »(°)] s! { [mpf {mr ] [t^wf {^ >(0) }1 } ' (3 - 36) 



where the set of partitions p(s, m) now contains only a single element 

p(s, m) = {(ci2) 04) : 02, «4 S A/o, a 2 + a 4 = m, 2a 2 + 4a 4 = s} . (3.37) 

Hence 

a-i = 2m — — s (3.38) 
1 

124 = — s — m, (3.39) 
for 4m > s > 2m and s even. Expression 13 . 36| is very simple to encode in a numerical subroutine. 

A. Boundary conditions 

Having separated Dirac equations in Anti-de Sitter Schwarzschild black hole background spaces and obtained two 
radial second order differential equations . 12 . 1 91 and 12.201 we have argued that those two equations are related. They 
can be physically identified as the equations governing the axial and polar perturbations of the Dirac field, respectively. 
We have also mentioned the fact that the two potentials 12 . 2 II and 12 . 221 possess the same quasinormal mode frequency 
spectra if the boundary conditions are mapped consistently. 

According to Cooper [23^. equation 12.191 describes axial perturbations and equation 12.201 describes polar pertur- 
bations of the Dirac field. Suppose we are interested in investigating further the AdS/CFT conjecture and want to 
relate the spinor field and its quasinormal frequencies to the calculation of the correlation functions on the boundary 
conformal field theory More precisely, let us suppose that we are interested in evaluating the correlators of a polar 
quantity in the CFT on the boundary, which corresponds to the polar part of the spinor field in the bulk of AdS 
space-time. Accordingly we must set the axial spinor field perturbation to vanish on the boundary. 

In terms of what we have done until now, we must take our series solution 13.221 to the axial equation 12.191 and 
impose on this Dirichlct boundary conditions 

00 

u{z) = a nz n -> (3.40) 

n=0 

as z — ► 1. Taking the limit of the left hand side gives us the identity 

00 

J2a n = 0. (3.41) 

n=0 

This is an implicit equation for the computation of quasinormal frequencies related to the spinor field. The computer 
program which performs the numerical calculation is organized in the following way. An independent routine evaluates 
the coefficients A s associated to the Taylor expansion of q~ 1//2 . These results are then exploited within a recursive 
routine which iterates the recurrence relation to evaluate the coefficients a n of the series solution. Finally, an ordinary 
Newton- Raphson root-finding routine seeks for solutions to the eouation l3.41l 

The choice of boundary conditions is related to the axial or polar character of the spinor field at infinity and to its 
interpretation in terms of the AdS/CFT conjecture. On the other hand the consistency of the boundary conditions 
can be checked by the shape of the effective potential in the differential eouation l2.19l we are solving. 

In Fig. ^ we plot the potential related to the evolution of the axial component of massless spin- 5 fields alongside 
the axial potential for the spin-2 case in four dimensional Schwarzschild- AdS space-time. Note that r — > 00 is a 
regular point of the differential eouation l2.19l and Dirichlet boundary conditions are a consistent choice. 



IV. NUMERICAL RESULTS 



The numerical procedure described above can be successfully applied to both large black holes, for which the radius 
of the event horizon n is much larger than the Anti-de Sitter radius a -1 (7*1 >> a -1 ), and intermediate black holes 
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u y O 2 4 6 8 10, 12 14 16 18 20 

FIG. 1: In this figure we show the potentials related to axial perturbations of AdS Schwarzschild black holes. The plot on the 
left represents the potential related to the evolution of an axial spin-5 field, for / = and M = 1.0. The plot on the right 
represents the potential related to the evolution of an axial spin-2 field, for I — 3 and M = 1.0. The masses are rescaled with 
respect to the radial coordinate of the black hole event horizons. 

for which r\ ~ a . The quasinormal frequencies are decomposed into real and imaginary parts 

u> = uji + iluj. (4-1) 
With this sign chosen, ioj is negative for all quasinormal frequencies. In Table [fl we list the values of the lowest 



ri 




LUI 


n 


LOR 


(j 1 


100 


~ 


-76.8157 


1.0 


1.80808 


-1.10565 


60 


~ 


-46.8901 


0.8 


1.76332 


-0.8516 


40 


~ 


-31.9303 


0.5 


1.7463 


-0.4525 


20 


~ 


-17.0652 


0.45 


1.7555 


-0.3829 


10 


~ 


-10.0749 


0.4 


1.7719 


-0.3118 


5 


2.20 


-6.3626 


0.3 


1.7958 


-0.1768 



TABLE I: Values of the lowest Dirac quasinormal mode frequency for a 1 = 1.0 and n = 0, 1 = 0. These results are obtained 
by evaluating the frequency as a function of r\ for some selected black hole sizes. 

quasinormal mode frequencies for aT x = 1, I = and selected values of 7*1. For large black holes, the real and 
imaginary parts of the frequency scale linearly with the horizon radius, resembling the results obtained in for 
other perturbation fields. Since the temperature scales also with r% in this regime, the imaginary part of the frequency, 
which determines how damped the mode is and which according to the AdS / C FT correspondence is a measure of the 
characteristic time r = l/|w/| of approach to thermal equilibrium, scales with the temperature. Therefore, in the dual 
CFT the approach to thermal equilibrium is faster for higher temperatures. The linearity of the scaling between w/ 
and T is clearly shown in Fig. |2 both for the n = and n = 1 mode numbers. The dots represent some frequencies 
numerically calculated. The lines connecting them are linear fits. Explicitly, the lines are given by 



-3.11T 



(4.2) 



for n = and 



-9.36T 



(4.3) 



for n = 1. 

In Table ^ we list the values of the first eleven quasinormal mode frequencies corresponding to massless spinor 
perturbations of large Schwarzschild-AdS black holes (ri ^> a" 1 ). For large values of the overtone number n, the 
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-90 | 

-80 - 
-70 - 

-60 - 

-40 - 

-30 - 

-20 - 

-10 - fi>'' 



FIG. 2: The imaginary part of the spinor quasinormal mode frequencies for n = 0, / = (left) and n = 1, I = (right) are 
shown as functions of the temperature T for some selected black hole sizes. The dots represent some frequencies numerically 
evaluated. The lines connecting them are linear fits. 



Overtone Number 


I -- 


= 


I -- 


= 1 


n 




(J I 


UJR 


UJI 





~ 


-76.8402 


~ 


-78.7731 


1 


~ 


-217.83 


~ 


-211.417 


2 


47.7354 


-417.775 


76.7222 


-406.192 


3 


178.073 


-650.494 


203.247 


-636.479 


4 


305.383 


-877.767 


330.457 


-863.679 


5 


433.279 


-1104.39 


458.292 


-1090.25 


6 


561.539 


-1330.64 


586.517 


-1316.47 


7 


690.04 


-1556.66 


715.01 


-1542.46 


8 


818.77 


-1782.56 


843.695 


-1768.29 


9 


947.614 


-2008.24 


972.531 


-1994.02 


10 


1076.56 


-2233.86 


1101.48 


-2219.65 



TABLE II: Quasinormal mode frequencies corresponding to I = and / = 1 massless Dirac perturbations of a large 
Schwarzschild-AdS black hole (ri = 100). For large values of the overtone number n, the modes become evenly spaced in 
mode number and the spacing is given by * — ~ (1.299 — 2.25i). 



-150 - 

3 

-100 - 



-50 - 



frequencies become evenly spaced in mode number and the spacing, which is independent of the angular mode number 
J, is given by 

(o; n+ i - u n ) ^ 9 _ 2 25 ., (4 4) 

r\ 

This is exactly the same spacing, related to the large black hole regime, obtained by Cardoso and his collaborators in 
|24| . for different kinds of perturbing fields. Their results were the same for scalar, electromagnetic and gravitational 
perturbations. Moreover, the quasinormal frequencies of large black holes have a number of first overtones with pure 
imaginary parts as in the electromagnetic and gravitational cases. The higher the black hole radius r\, the higher the 
number of these first pure damped modes. In Fig. [21 the large black hole frequencies are plotted for I = 0, 1. 

In Table ITTT1 we list the first ten quasinormal mode frequencies corresponding to massless spinor perturbations of 
intermediate Schwarzschild-AdS black holes (n ~ a -1 ). The values of the frequencies are shown for a~ l = 1.0 and 
I = 0, 1. For large values of the overtone number n, the modes become evenly spaced and the spacing, which is 
independent of I , is given by 

n 

Once again this is a result that, for the intermediate black hole regime, resembles those obtained in [24| for different 
kinds of perturbation fields. In Fig. 0] the intermediate black hole frequencies are plotted for I = 0,1. 
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J 
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UJR 

FIG. 3: The frequencies of the first eleven quasinormal modes for the axial component of the Dirac field, are shown for a" 1 = 1, 
ri = 100 and the angular mode values I = and 1 = 1. 



Overtone Number 


I 


= 


I -- 


= 1 


n 




m 


LOR 


UJI 





1.80808 


-1.10565 


2.7956 


-0.9744 


1 


3.44229 


-3.43805 


4.26976 


-3.1623 


2 


5.28094 


-5.82213 


5.99062 


-5.51618 


3 


7.17816 


-8.119914 


7.81556 


-7.89152 


4 


9.09828 


-10.5728 


9.68897 


-10.2668 


5 


11.0301 


-12.9413 


11.5837 


-12.6412 


6 


12.9712 


-15.3067 


13.5014 


-15.0117 


7 


14.9172 


-17.6699 


15.4250 


-17.3815 


8 


16.8671 


-20.0254 


17.3552 


-19.7468 


9 


18.8253 


-22.3793 


19.2938 


-22.1086 



TABLE III: Quasinormal mode frequencies corresponding to I = and I = 1 massless Dirac perturbations of an intermediate 
Schwarzschild-AdS black hole (ri = 1.0). For large values of the overtone number n, the modes become evenly spaced in mode 

number and the spacing is given by '""+,! ~"" ) ~ (1.96 - 2.36i). 

In Fig. and we show, for intermediate black holes, the dependence of the real and imaginary part of the 
fundamental quasinormal mode frequencies on the metric parameters for I = and 1 = 1. The parameters k and a" 1 
represent, respectively, the black hole surface gravity evaluated at the event horizon and the Anti-de Sitter radius. 
Both the real and imaginary parts of the frequencies seem to depend linearly on the black hole surface gravity and 
therefore on the black hole temperature for > 1.8. This is a result that has already been found in [7j for 

gravitational perturbations of Schwarzschild AdS black holes. The slopes of the straight lines, which are independent 
of /, are given by 



-2000 
-1500 

3 

-1000 
-500 



lur = 0.25 k 



(4.6) 
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FIG. 4: The frequencies of the first ten quasinormal modes for the axial component of the Dirac field, are shown for n = a -1 = 
1.0, and the angular mode values I = and I = 1. 

for the plots in Fig. [5] and 

u! T = -1.29k (4.7) 

for the plots in Fig. El 

For small black holes (ri <C ot ), it is very difficult to evaluate the corresponding quasinormal frequencies, or to go 
high in mode number. The error associated in estimating the frequencies in this regime is too high, and we cannot be 
completely sure of the results. In the case of other perturbation fields 0,E|I3]j ^ was found that small black holes 
have quasinormal frequencies that are very close to the pure AdS values. In the case of the Dirac field considered 
in this paper, after setting a -1 = 1, we have been able to evaluate the fundamental frequency, at most, down to the 
value of r*i = 0.3 (see Table QJ. Further investigations and more powerful numerical tools are needed to outline the 
behaviour of the spin-1/2 frequencies in the limit of r% — > 0. 

V. CONCLUSION 

In this paper we have investigated the quasinormal modes frequencies of a spinor field interacting with a 
Schwarzschild black hole in four dimensional Anti-de Sitter space-time. After solving the related axial perturbative 
equation we have imposed on this solution Dirichlct boundary conditions and evaluated the corresponding quasinormal 
frequencies. This choice of boundary conditions is physically justified in the contest of the AdS/ CFT correspondence, 
in which it is suggested to exploit these frequencies in the evaluation of poles of correlation functions associated to 
the conformal field theory on the boundary. 

For large black holes, we have numerically confirmed that the imaginary part of the frequency scales with the 
temperature. We have also shown that both the real and imaginary components of the frequencies are evenly spaced 
in mode number. The spacing between consecutive modes, which is independent of the angular quantum number I, 
behaves as in the case of scalar, electromagnetic and gravitational perturbations. The quasinormal frequencies have 
a number of first overtones with pure imaginary parts. 

For intermediate black holes, we have found that both the real and imaginary part of the frequency scale linearly 
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FIG. 5: The dependence of the real part of the fundamental quasinormal mode frequencies on the metric parameters for / = 
(bottom) and 1=1 (top), is shown for Dirichlet axial boundary conditions on the spinor field. 



with the surface gravity. The modes are evenly spaced and the spacing, which is independent of the angular quantum 
number I, resembles the one related to other kinds of perturbations. 

Our purpose is to extend this project imposing different boundary conditions on the axial or the polar Dirac 
equations. We could be interested, for instance, in calculating the correlation functions and associated poles of a 
quantity related to the axial part of the spinor field. That would suggest we set the polar component of the spinor 
field perturbation to vanish on the boundary and consequently we could solve either the odd parity differential 
equation with Dirichlet boundary conditions or the even parity equation with mixed boundary conditions due to the 
transformation theory between polar and axial functions. 
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